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Since the very first experiments, superconducting circuits have suffered from strong coupling to 
environmental noise, destroying quantum coherence and degrading performance. In state-of-the-art 
experiments, it is found that the relaxation time of superconducting qubits fluctuates as a function 
of time. We present measurements of such fluctuations in a 3D-transmon circuit and develop a 
qualitative model based on interactions within a bath of background two-level systems (TLS) which 
emerge from defects in the device material. In our model, the time-dependent noise density acting on 
the qubit emerges from its near-resonant coupling to high-frequency TLS which experience energy 
fluctuations due to their interaction with thermally fluctuating TLS at low frequencies. We support 
the model by providing experimental evidence of such energy fluctuations observed in a single TLS 
in a phase qubit circuit. 
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I. INTRODUCTION 

Superconducting qubits 1 are well on the way towards 
achieving the prerequisites for fault-tolerant quantum 
computation schemes 2-4 . With the advent of highly 
coherent superconducting circuits for quantum applica¬ 
tions, previously neglected sources of environmental noise 
become important. One such cause of decoherence is 
spurious two-level systems (TLS), which are believed to 
be present in large numbers in the amorphous dielectric 
oxide layer covering the superconducting films 5,6 . En¬ 
sembles of TLS naturally explain the low-temperature 
properties of glasses 7,8 and are used as a universal model 
for l//-type low-frequency noise in electric circuits 9,10 . 

Virtually all designs of superconducting qubits tested 
so far show a pronounced frequency dependence in their 
relaxation rates 11-14 , which indicates strongly coloured 
high-frequency noise acting on the circuits 15 . A natural 
explanation of these observations relies on weak inter¬ 
actions between the circuit dynamics and spurious en¬ 
vironmental TLS, possibly located in the disordered di¬ 
electric covering these circuits. For coupling strengths 
that are much weaker than the individual decoherence 
rates of qubit and defect, the effect of the TLS on the 
qubit will be that of a strongly peaked high-frequency 
noise spectrum. In other experiments, strongly coupled 
coherent TLS are often found to cause avoided level cross¬ 
ings in superconducting circuits which include Josephson 
junctions 5,6 . Those TLS are believed to reside in the di¬ 
electric forming the tunnelling barrier inside the circuits 
Josephson junctions, enabling their stronger coupling to 


the circuit dynamics. Otherwise they are conjectured to 
be of the same origin as the TLS observed as resonances 
in the high-frequency noise spectrum. Using supercon¬ 
ducting qubits as probes, it is possible to fully charac¬ 
terise the properties of the strongly coupled defects, for 
example by measuring their level-structure and coherence 
times 16-18 . 



Figure 1. (Color online) Fluctuations in time of the relaxation 
rate Ti of a superconducting 3D-transmon qubit. Errorbars 
show the 95% confidence interval of the fits, the dotted red 
line indicates the mean value of all measurements and the 
dashed black line is a moving average over 10 samples to em¬ 
phasize the multivalued character of the jumps. Each individ¬ 
ual point in this plot required a measurement time of ~ 1 min. 


In this work, we discuss the origin of time-dependent 
fluctuations in the energy relaxation time T), which 
are observed in superconducting 2D-transmons 19 , flux 
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qubits 20 and 3D-transmons 14 , as shown in Fig. 1. The 
paper is organised as follows: Section II starts by mo¬ 
tivating this work and our theoretical model for time- 
dependent fluctuations in the relaxation rate of super¬ 
conducting circuits. Section III then describes the ex¬ 
periments from which our data originate. In the follow¬ 
ing part, Sec. IV, we develop the model and present the 
main results. The discussion in Sec. V presents implica¬ 
tions and possible tests of the model and considers pos¬ 
sible alternative explanations of the data. The paper is 
followed by appendices summarising details of the experi¬ 
ments, additional experimental data, and providing more 
details on the theoretical calculations. 


II. MOTIVATION 

Qubit relaxation may occur through its weak coupling 
to environmental TLS whose characteristic eigenenergies 
are comparable to the qubit’s energy splitting. The en¬ 
vironmental noise spectral density originating from cou¬ 
pling to a single such TLS is strongly peaked around its 
eigenfrequency. A natural approach to explain the fluc¬ 
tuations in the qubit relaxation rate is thus to assume 
random changes in the energy splitting of individual two- 
level defects; c.f. Fig. 2. Our model for the origin of 
the fluctuations is then based on the presence of a large 
number of interacting TLS at both low and high eigenfre- 
quencies. Due to the interactions between TLS, thermal 
switching of the state of low-frequency TLS will then lead 
to fluctuations in the energy splitting of high-frequency 
TLS, providing a qualitative description of the observed 
data. This model is further underpinned by our direct ob¬ 
servation of fluctuations in a high-frequency TLS’ energy 
splitting, which occurs on time scales comparable with 
the qubit’s Tj fluctuation; see Fig. 3. In the following, we 
will indicate TLS with eigenenergies much larger than the 
thermal energy as TS (tunnelling systems), while those 
at energies much lower than temperature will be named 
TF (thermal fluctuators). 

Our model provides a qualitative description of the 
origin of fluctuations in the electrical susceptibility of 
mesoscopic circuits, an area which has recently started to 
attract attention from both experiment and theory 21-23 . 
We also note that interactions between TLS have recently 
been observed directly in two strongly coupled defects 24 
and that such coupling has been invoked as a model of 
noise before, e.g. to explain the line-width broadening 
and spectral diffusion of ultrasonic excitations of TLS en¬ 
sembles in glasses 25,26 as well as spectral blinking of dye 
molecules 27 and quantum dots 28 . More recently, Refs. 29 
and 30 make a connection between slow fluctuations in 
the resonance frequency of superconducting resonators 
causing phase noise, and ensembles of interacting TLS 
leading to fluctuations in the energy splitting of high- 
frequency TS, much along the same lines as we describe 
here. While in that work the real part of the suscepti¬ 
bility was considered, leading to fluctuations in the level 



Figure 2. (Color online) Illustration of the conjectured mech¬ 
anism behind fluctuations in IT. We plot the noise spectral 
density C(u), Eq. (2), of a single high-frequency TS as a func¬ 
tion of frequency u. The qubit level splitting is indicated as 
uj q and the TS energy as E. Fluctuations in E, as indicated by 
the arrow and the dashed contours, may cause strong changes 
in the noise spectral density at the qubit frequency, leading to 
significant changes in the qubit relaxation rate Fi oc C(oj q ). 
The inset shows an illustration of the interaction between a 
central high-frequency TS (red, centre) with a surrounding 
bath of low-frequency TF (black), where the interaction is 
limited to a small spatial range indicated by the grey-shaded 
region. 


splitting of a resonator, here we are concerned with its 
imaginary part that is responsible for energy dissipation. 

III. EXPERIMENTAL EVIDENCE 

The fluctuations of the Tj-time reported here (Fig. 1) 
were measured in a superconducting qubit in the 3D- 
transmon design 14 , with an average relaxation time T) 
of ~ 80 /zs. In our 3D-transmon circuit, the qubit en¬ 
ergy, i.e. the level splitting of its two lowest levels, is 
fixed at ui q /2 tt = 3.58 GHz and not tuneable as in other 
designs 31-35 . Each datapoint results from a series of indi¬ 
vidual measurements, each time resonantly exciting the 
qubit and detecting the qubit population after waiting 
for some time t. The resulting traces where fitted to an 
exponential decay curve oc e -rit . The observed fluctua¬ 
tions of the qubit’s relaxation rate Fi do not show any 
apparent structure, with the largest experimentally re¬ 
solvable fluctuation rate given by the inverse of the time 
it takes to obtain a single value of T), here ~ 1 min. 
Additional datasets are shown in Appendix B. 

In a second experiment, we use a superconducting 
phase qubit to directly monitor the properties of a single 
high-frequency TS that is strongly coupled to a super¬ 
conducting phase qubit. Figure 3 (a) shows measured 
time-dependent fluctuations of the TS’ energy level split¬ 
ting which occur on time scales similar to those of the 
above discussed qubit fluctuations. Here, the TS’ res¬ 
onance frequency E was repeatedly measured by vary- 
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Figure 3. (Color online) (a) Change in TS energy E as a 
function of time measured in a superconducting phase qubit. 
Errorbars indicate the 95% confidence interval of the fits, the 
red dotted line is the average over the samples shown and 
the black dashed line is a moving average over 10 samples. 
Panel (b) shows two Lorentzians in the escape probability of 
the qubit P(|l)) at two different times as an example of the 
change in TS energy. Here the dots are the raw data and the 
solid lines are the result of a fit to the data. Vertical dashed 
lines in (a) indicate the measurement times for the two curves 
shown. 


ing the frequency of a long microwave pulse applied to 
the qubit circuit with a pulse amplitude that was large 
enough to allow for direct excitation of the TS excited 
state 1 '. During the microwave pulse, the qubit was kept 
far detuned from the TS. After the pulse, qubit and TS 
were brought into resonance in order to swap the TS ex¬ 
citation into the qubit, whose population was then read 
out. Details of this technique can be found in Refs. 17 
and 24 and Appendix A. 


IV. THEORETICAL MODEL 
A. TLS as sources of fluctuating noise 


In the following, we describe our model explaining the 
observed fluctuations in the relaxation rate Ti = 1/T) of 
superconducting circuits. We first note that in a mas¬ 
ter equation description of dissipative quantum dynam¬ 
ics, the relaxation rate of a qubit is proportional to the 
unsymmetrized spectrum of its environment at the fre¬ 
quency of the qubit’s level-splitting, Ti cx C(w 9 ) 36 . Here 


we assume effectively zero temperature, ksT -C hui q , so 
that thermal excitations can be neglected. It is then our 
goal to relate fluctuations in the energy of a single TS to 
changes in the high-frequency noise acting on the super¬ 
conducting circuit and to further characterise the fluc¬ 
tuations in terms of parameters of the experiments and 
the TLS distributions. We start by describing a single 
TLS as a quantum two-level system using the tunnelling 
Hamiltonian 7,8 

#tls = - \ ea * + > W 

where e is the asymmetry energy between the two wells 
and A is the tunnel splitting. The Pauli-matrix <j z here 
describes the position of a particle on either side of a 
double-well potential, and a x induces transitions between 
the wells. Diagonalizing yields Htls = \E& Z with the 
level-splitting E = y/ e 2 + A 2 . Here and in the following 
we use the convention h = 1 , so that all energies are 
expressed in units of angular frequencies. 

The TS observed in high-frequency noise spectra are 
believed to be charged entities interacting with the su¬ 
perconducting circuits via their electric dipole moment 
oc cr 2 5 . Assuming weak qubit-TS coupling, their effect on 
the qubit will be that of strongly coloured noise, where 
the spectral density can be calculated from the Fourier 
transform of the two-time correlation function of their 
coupling operator a z 15 . We obtain 


C(w) = J dte lut (cr z (t)cr z (0)) 
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with the TLS’ equilibrium steady-state population 
(a z ) = tanh (E/2J:bT), the intrinsic TLS relaxation rate 
71 , and 72 = §71 + 7 ^,, where 7 ^ is the pure dephas¬ 
ing rate of the TLS. Here, tanf? = A/e defines the TLS’ 
mixing angle. Eq. (2) is composed of three parts, each of 
which is relevant for TLS in different parameter regimes. 
The first line describes low-frequency noise due to ran¬ 
dom switching of the TLS and is most pertinent for low- 
frequency TF with E < ksT. The second term is a high- 
frequency contribution which is sharply peaked around 
the TLS energy and is most pronounced for TS with 
E ksT. Since those TS are mostly resting in their 
ground state, they are able to absorb energy from the 
qubit. It is this contribution that gives rise to the ob¬ 
served resonances in the noise spectrum 11-13 and in which 
we are mostly interested. The final term contributes at 
negative frequencies and describes the ability of the TLS 
to excite the qubit by transferring an excitation to it. For 
both high-frequency TS in thermal equilibrium as well as 
low-frequency TF, this term will not contribute. 
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For simplicity, we assume the environmental noise at 
frequencies close to the qubit level splitting to q is domi¬ 
nated by a single, weakly coupled high-frequency TS at 
energy E ~ u> q . We further assume this TS is interacting 
with a large number of other TLS which are located in 
its close spatial vicinity. This is the situation illustrated 
in Fig. 2 and the one most relevant to experiment 11-13 . 
If the distribution of TS at high frequencies is dense 15,29 , 
our results still hold but have to be additionally averaged 
over the high-frequency distribution. We model the inter¬ 
action between all TLS in the sample by a Hamiltonian 
of the form 

( 3 ) 
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where g :/ is the coupling strength between the high- 
frequency TS and all other TLS, indicated by the index 
j. Coupling of the type Eq. (3) can be caused e.g. by 
electric dipole coupling or strain-mediated interaction, 
where the asymmetry-energy of either TLS depends on 
the relative position of the other TLS in their respective 
double-well potentials 24 . With such an interaction, the 
energy splitting of any TLS depends on the instantaneous 
state of all TLS in a certain range around it, determined 
by the microscopic origin of their interaction; c.f. inset 
to Fig. 2. 

We are looking at fluctuations in the qubit relax¬ 
ation rate due to slow fluctuations in the TS energies 
E. In order to calculate expectation values and statis¬ 
tics, we write the level splitting of an individual TLS in 
the form E = E 0 — . gja z j, now depending on the 

state of all other TLS via the mutual interaction g-j from 
Eq. (3). Here we focus on a high-frequency TS with 
E 0 ksT, 72 , such that (a z ) = — 1 and the result¬ 
ing spectral density is strongly peaked around the TS 
eigenenergy E, c.f. Eq (2). We defined the undisturbed 
TS level splitting as Eq = y/£ 2 + A 2 with the parameters 
e and A from Eq. (1). 

We can further write the qubit relaxation rate due to 
its coupling to a single high-frequency TS as Ti oc 7 g . 
Here, the relaxation rate induced by a single TS is given 
by the high-frequency components of its spectral density, 
c.f. Eq. (2), as 
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Assuming the interaction between individual TLS to be 
weak, gj <C 72 , we can expand this to first order as 
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Equations. (5)-(7) will be the basis for our further calcu¬ 
lations. 


B. Distribution of TLS parameters 

For tunneling TLS one usually assumes flat distribu¬ 
tions for both the asymmetry energy e as well as the 
tunneling barrier height 7,8 . Since the tunneling energy A 
depends exponentially on the barrier, the resulting dis¬ 
tribution in TLS parameters is P(e, A) ~ 1/A. The TLS 
relaxation rates are then also distributed log-uniformly, 
P( 71 ) ~ I/ 71 , since the tunnelling strength depends 
mainly on the size of the tunnelling barrier. In Ref. 15 
it was found that a linear or super-linear distribution in 
e would naturally explain both low- and high-frequency 
parts of the noise spectrum acting on the qubit as stem¬ 
ming from the same ensemble of TLS. For the sake of 
generality, we will therefore assume the distribution of 
TLS parameters as 


P(e, A)d£dA = A ^-cfedA , ( 8 ) 

with a > 0 and the constant A needed for normalization. 
For non-interacting TLS, the distribution is usually as¬ 
sumed to be flat, a = 0 7,8 , but might be different from 
zero in the more realistic case of interacting TLS 15,29 . 
Without loss of generality we restrict the integration to 
the positive real axis. The distribution of inter-TLS cou¬ 
pling strengths gj depends strongly on the physical model 
of their interaction. It is important to note that the cou¬ 
pling strength g in most models can be both positive or 
negative, meaning the coupling between the TLS can ei¬ 
ther raise or lower the energy of the respective partners. 
For the dipole coupling model this reflects the fact that 
the relative orientation of the dipoles can be both parallel 
as well as antiparallel. 


C. General considerations 

In the calculations one has to carefully separate the 
different time scales of the problem. The measurement 
protocol fixes three distinct scales, which have to be com¬ 
pared to the fluctuation rates of individual low-frequency 
TF to determine the nature of their contribution to the 
fluctuations in the qubit’s relaxation rate F|. First, there 
is the time it takes to do a single measurement of the 
qubit population, t meas , where many such measurements 
are averaged to obtain each point in a complete relax¬ 
ation curve. Fluctuating TF that are faster than l/t meas 
will not contribute since they average out even for a sin¬ 
gle measurement. Second, there is the time to measure 
a single point of a curve, f po int- Fluctuations that are 
faster than l/tpoint, but slower than l/t meas w iH act as 
an effective broadening of the high-frequency TS reso¬ 
nance, increasing its line-width 72 . The slowest time scale 



5 


is given by the duration of the measurement of a com¬ 
plete Tj curve, tr, ■ TF dynamics slower than l/f po int 
but faster than l/tr 1 will lead to jitter in the energy re¬ 
laxation curve, contributing additional noise in the fit of 
Tj. Finally, slow TF that fluctuate at frequencies that 
are smaller than 1/t t x will be the ones responsible for 
the low-frequency fluctuations visible in the Xj data, see 
Fig. 1. Note that the microscopic origin of these small 
switching rates is so far unclear 29 . For the very slow fluc¬ 
tuations observed in experiments, on time scales ~ min, 
to the best of our knowledge no microscopic model ex¬ 
ists. A possible candidate might be collective behaviour 
of large ensembles of TLS that form clusters 3 '’ 38 , but 
clear experimental confirmation of this effect is missing 
so far. 

In the following we will be interested in calculating 
the temperature and frequency dependence of the qubit 
relaxation rate due to its coupling to individual TS, av¬ 
eraged over TLS parameter distributions, as well as the 
spectrum of the fluctuations in Ti. Due to the considera¬ 
tions above, the temperature dependence will be strongly 
influenced by the thermally activated part of the TF dis¬ 
tribution, contributing via the TS linewidth 72 . Follow¬ 
ing Refs. 15 and 29 we find the temperature dependence 
of the dephasing rate due to a bath of low-frequency TFs 
as 72 oc T a+l , where a characterizes the TLS distribu¬ 
tion. 


D. Average relaxation rate 

We now turn to calculating the average of the qubit’s 
relaxation rate using the distributions introduced above. 
We concentrate here on fluctuations originating from the 
low frequency contributions from TF with small level 
splitting, E < UbT, since those are the ones directly ob¬ 
servable in experiment. Noting that for TLS in thermal 
equilibrium (d z ) = cos6(a z ) = cosdtanh (E/2IcbT), we 
can directly write down the mean value of the qubit relax¬ 
ation rate due to the high-frequency TS to lowest order 
in the inter-TLS coupling strength g as 

(7?) = 7 q 0) + 7g 1} g i cos 6 J tanh ( 9 ) 

3 

where the sum includes all other two-level defects that a 
single high-frequency TS interacts with. In the calcula¬ 
tion of the average rate ( i y q ), we immediately notice that 
f dg gP(g) = 0 , since we integrate an odd function over 
an even range. Therefore we simply find 

(%) = 7g 0) j (10) 

i.e., the average contribution to the relaxation rate from 
a single TS is given by it’s spectrum centred around its 
undisturbed level-splitting E 0 . For the temperature de¬ 
pendence of the ensemble-averaged qubit relaxation rate 


we then find 

/T ^ 272 

< r ’> “ 


6 uj < 72 
Suj » 72 


( 11 ) 


where Suj = uo q — E is the detuning between qubit and 
TS, and we distinguish between the case where qubit and 
TS are nearly resonant and when they are far detuned. 


E. Rate fluctuation spectrum 

The spectrum of fluctuations of the qubit relaxation 
rate is then related to the Fourier transform of the rate 
correlation function as 

(7g( i )7g(°)) w = f dt (% W^ 0 )) 

= ( 7 9 1} ) J29j9l(^z,j(t)a z ,i{0)) u ( 12 ) 

3,1 

where, in evaluating the correlator, we restrict ourselves 
to the low frequency contribution of the TLS autocor¬ 
relation function Eq. (2), i.e. we focus on TFs with 
E <C ksT. Additionally we assume that different TLS 
are uncorrelated, (a z ja Zt i) = 0. We are also only inter¬ 
ested in the bare fluctuations of the rate, so we have al¬ 
ready subtracted the mean value above. For more details 
on the calculations, see Appendix C, where we addition¬ 
ally discuss the case when the TF switching is solely due 
to interactions with phonons. 

For the average over the coupling strength, one finds 
J dg g 2 P(g) oc const, where the constant is mainly deter¬ 
mined by the maximum possible coupling strength and 
thus by the minimal distance between TLS and the mi¬ 
croscopic origin of their interaction. Performing the aver¬ 
age over the mixing angle 6 also contributes a constant, 
with the exact value again depending on details of the 
microscopic TLS model. The average over TF energies 
can be written as 

/ dE p < £ > 0 - ta " h2 (^f))*£ iEEa = T ° + ' 

(13) 

contributing to the temperature dependence of the final 
result. Still assuming small interaction strength between 
TLS, g <C 72 , we can now distinguish three regimes re¬ 
lated to the initial detuning between our qubit and the 
high-frequency TS, Suj = u> q — Eq- For qubit and high- 
frequency TLS near resonance, Suj <C 72 , we find that 
7 q 1 ^ oc Suj/ 7 I, while in the regime of intermediate detun¬ 
ing, Suj ~ 72 , one finds 7 g 1 ' ) oc I/ 7 I. In the far detuned 
regime, Suj ^$> 72 , we finally have jij 1 ' 1 oc 72 /Suj 3 . Finally, 
adopting the standard assumption for tunnelling TLS, 
P( 7 i) ~ I/ 71 , the frequency dependence of the fluctua¬ 
tion spectrum is determined by 

/■™ ax 27 ! 1 

/ <^7l P{ll)~ 2 ~, -7 OC - , UJ < 7Max ■ (14) 

Jo 7f + w- UJ 
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Here the maximum relevant switching rate 7 Max is given 
by the time of a single Tj-measurement. All faster fluc¬ 
tuations will be averaged out in the observations, leading 
to the behaviour ~ 1 /w for ui < 7 Max- 

Thus, we find the temperature and frequency depen¬ 
dence of the Ti fluctuation-spectrum as 

f T - 5(a+i) 5 Su < 72 

(r 1 (t)r 1 (0)) u cxu J - 1 \T- 3 ^ , 6u~ 72 , (15) 

[t 2 (“ +1 ) , 5u> > 72 

where (.. .) w denotes the Fourier transform of the two- 
time correlation function of the relaxation rate T \. These 
results hold for small inter-TLS coupling g :l <C 72 - The 
opposite case gj 72 corresponds to on-off switching 
and is excluded by the experimental data. In the inter¬ 
mediate regime, gj ~ 72 , the overall temperature depen¬ 
dence will be given as an average over our results. 


V. DISCUSSION 

A. Implications and tests of the model 

Our model can be directly tested by measuring the 
relaxation rate at different qubit level splittings and in¬ 
ferring the time and frequency dependence of the noise 
spectrum acting on the qubit. By using a frequency- 
tunable qubit, the fluctuations in the noise spectral den¬ 
sity might be directly resolvable in time and frequency, 
depending on the time-scale of a single measurement of 
the relaxation time Tj. Even for non-tuneable qubits it 
is possible to probe the noise spectral density in close 
vicinity of the qubit frequency by measuring the decay 
of Rabi oscillations of the qubit, c.f. Appendix D and 
Ref. 39. Another possibility is to apply external driving 
to saturate the TF responsible for the fluctuations in TS 
energy. If an electric field is applied resonantly with the 
relevant low-frequency TF, it will lead to oscillations with 
the Rabi frequency depending on the detuning between 
drive tone and TF energies, the TF dipole moments and 
the electric field strength at their position. Assuming 
the resulting Rabi frequency is fast compared to the du¬ 
ration of a single X| measurement, the effect would be 
to raise the average (Ti) while at the same time reduc¬ 
ing the amplitude of its fluctuations. This is because 
the resonant driving of initially very slow TF will alter 
their contribution towards a simple line-width broaden¬ 
ing of the high-frequency TS. In experiments with 3D- 
transmon qubits this could be achieved by careful engi¬ 
neering of the cavity modes, such that there exists a suit¬ 
able low-frequency mode exhibiting strong electric field 
components spatially close to the qubit. In other qubit 
architectures this might be possible within existing ex¬ 
perimental setups 17 . In our transmon qubit sample, this 
experiment proved unfeasible due to design restrictions 
in the employed cavity. Additional verification could be 


achieved by a systematic characterisation of the fluctua¬ 
tions of T\ at a variety of experimental temperatures T. 
An additional challenge arrises from the fact that the ex¬ 
act temperature dependence is connected sensitively to 
the qubit-TS detuning Soj, c.f. Eq. (15), which also has 
to be determined in this case. 


B. Alternative models 

Possible alternative models for the fluctuating noise 
spectrum include fluctuations of the quasiparticle density 
in the superconductor. Quasiparticle tunnelling across 
the circuit’s Joseplrson junctions can induce relaxation 
and dephasing 40 , and explains well the temperature de¬ 
pendence of qubit relaxation rates for elevated sample 
temperatures. In contrast to our model, which depends 
on a structured noise spectrum as background, the quasi¬ 
particle induced noise is flat at high-frequencies. Follow¬ 
ing Ref. 40 we calculate the fluctuations in quasiparticle 
density required to effect the observed variance in the re¬ 
laxation time of transmon qubits. For the parameters of 
our sample, we find the fluctuation in the quasiparticle 
volume density required to change the relaxation rate by 
1 kHz as Sn qp « 0.5/^rn 3 , see Appendix E for details. 
From the geometry of our sample, it then follows that 
this change would require the number of quasiparticles 
present on either one of the qubit islands to fluctuate 
by SN qp « 1.5 x 10 4 . We are not aware of any mecha¬ 
nism leading to symmetric fluctuations in the quasipar¬ 
ticle number of this magnitude. 

Another possible model is that in the 3D-transmon 
sample used to obtain Fig. 1, the qubit level splitting 
might fluctuate in time, e.g. due to changes in the crit¬ 
ical current of the circuits Joseplrson junction 41,42 . To¬ 
gether with the observed strong structure in the noise 
spectrum 11-13 this could also explain the fluctuations in 
the qubit relaxation. Here we again have to be mind¬ 
ful of the time scales involved. Fast fluctuations of the 
qubit energy, i.e. faster than the Rabi frequency used to 
excite it (in our experiments tt/2n ~ 4 MHz), will not 
lead to the observed slow fluctuations in the relaxation 
rate, but rather average out over the measurement time 
fr, • Their effect would be such that the observed qubit 
relaxation rate would no longer depend on the noise spec¬ 
trum at a single frequency, but rather an average over the 
spectrum at a range of frequencies. Intermediate energy 
fluctuations, faster than l/t-Ti but slower than the Rabi 
frequency, will lead to a different qubit energy at each 
point of a single measurement and thus manifest as an 
additional source of noise in the fit parameters. Very 
slow fluctuations in the qubit energy, on the same time 
scales as the variations in Xj will, in addition to poten¬ 
tially impacting the qubit relaxation rate, also influence 
the resonance condition for the Rabi pulse used to excite 
the circuit for each measurement. From our experimental 
data we obtain the excitation amplitude as a fit parame¬ 
ter for each measurement, c.f. Appendix B. While there 
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are fluctuations visible in these parameters, they are gen¬ 
erally of small amplitude and not highly correlated with 
the observed T) variations. We therefore conclude that 
while this mechanism might be present, its effect is likely 
to be smaller than the one associated with fluctuating 
TLS. 


VI. CONCLUSION 

In this paper we present a simple model of interact¬ 
ing TLS which offers a qualitative understanding of the 
observed fluctuations in the relaxation times T\ of super¬ 
conducting quantum circuits. The model is grounded in 
our experimental observations, grants a clear route to¬ 
wards further confirmation, and provides a way to verify 
and refine the existing microscopic TLS models. More¬ 
over, our model clearly indicates that parasitic TLS are a 
limiting factor for the stability of today’s best performing 
superconducting circuits. A better understanding of this 
decoherence source is thus vital for further improving the 


fidelity of superconducting quantum circuits. 
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APPENDIX A: SAMPLE PARAMETERS AND 
EXPERIMENTAL REMARKS 

A. Observation of Ti-fluctuations of a 3D-transmon 

The three-dimensional cavity resonator used in this 
work was machined from bulk aluminium 6061. The cav¬ 
ity has a nominal size of 18.6 mm x 15.5 mm x 4.2 mm 
engineered to give a resonant frequency of approximately 
12 GHz. Two bulk head SMA connectors are used as 
input and output ports. The loaded quality factor of 
the waveguide cavity is 3 x 10 4 , with the output connec¬ 
tor stronger coupled than the input connector in order 
to guarantee a high signal-to-noise ratio. The sample is 
shielded with a cryoperm can that is thermally anchored 
to the mixing chamber of the dilution refrigerator. 

The qubit manipulation and readout pulses are deliv¬ 
ered to the cavity via a single coax line filtered by a 10 dB 
attenuator at each temperature stage of the refrigerator. 
Input and output ports are directly connected to SMA 
filters made of eccosorb in order to block infrared radia¬ 
tion and thermalize the central conductor of the coupling 
connectors. 

The output readout signal is amplified by a chain of 
cryogenic and room temperature amplifiers for a total 
gain of 60 dB. A low pass filter and two cryogenic cir¬ 
culators are used between the sample and the cryogenic 
amplifier. The qubit state is readout via the dispersive 
shift of the waveguide cavity 43 . 

The qubit is fabricated on sapphire substrate via alu¬ 
minium double angle evaporation. Two rectangular pads 
of 350 /iin x 700 /iin, separated by 50 /an, are connected 
by a 1 /iin wide aluminium strip with a Josephson junc¬ 
tion of size 0.1 /an x 0.1 /am. The chip has a total size of 
3.0 mm x 6.7 mm and is kept in place in the waveguide 
cavity by small pieces of indium. The qubit is placed at 
the maximum of the electric field of the first cavity mode. 

The qubit used in this work has an energy gap w 9 /27r of 
3.5825 GHz, and anharmonicity (~ Eq/ 27t) of 171 MHz. 
The qubit is designed to work in the transition regime, 
with Ej/Eq ~ 61 44 . 

For the Tj measurements reported here, the qubit was 
resonantly excited with a microwave pulse amplitude 
leading to the Rabi frequency of H/27T = 3.5714 MHz, 
corresponding to a 7r-pulse duration of 140 ns. 

B. Observation of TLS frequency fluctuations 

Our theoretical model to explain time-dependent fluc¬ 
tuations in the energy relaxation time of superconduct¬ 
ing qubits is based on their near-resonant coupling to 
high-frequency TLS, with the additional assumption that 
those TLS themselves experience resonance frequency 
variations due to their interaction with thermally fluc¬ 
tuating defects at low frequencies. In this work, we in¬ 
clude first experimental evidence that individual high- 
frequency TLS may indeed show resonance frequency 
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fluctuations in time as shown in Fig. 3. 

In order to access TLS individually, we exploit their 
strong coupling to the state of a superconducting phase 
qubit when they are residing in the amorphous tunnel 
barrier of the qubit’s Josephson junction. We were us¬ 
ing a phase qubit sample that has been developed in 
the group of Prof. J. Martinis at University of Cali¬ 
fornia, Santa Barbara, USA, with sample parameters as 
described in Ref. 45. 

We recorded the Lorentzian resonance curve of the 
TLS by varying the frequency of a long microwave pulse 
applied to the circuit while the qubit was kept far de¬ 
tuned from the TLS. As described in Ref. 17, this allows 
one to resonantly drive TLS while they remain effectively 
decoupled from the qubit dynamics. To read out the TLS 
quantum state, the qubit is first prepared in its ground 
state and then tuned into the TLS resonance. This re¬ 
alises an iSWAP operation that maps the TLS state onto 
the qubit, where it can be measured. 

Some of the TLS that were investigated with this 
method showed time-dependent fluctuations of their res¬ 
onance frequency that were large enough to be resolved 
spectroscopically. Often, we observe telegraph-signal like 
switching of TLS resonance frequencies between two sim¬ 
ilar values, indicating coupling to one dominating ther¬ 
mally activated TLS at low frequency. 

To characterise the internal TLS parameters, tun¬ 
nelling energy A and asymmetry energy e were measured 
by recording the strain dependence 46 of its resonance fre¬ 
quency and performing a hyperbolic fit to the equation 
E = \/A 2 + e 2 . Figure 3 was obtained on a TLS that 
had A/2-k = 7.056 GHz and whose asymmetry energy 
was tuned to e/27r = 918 MHz. At this asymmetry, this 
TLS had an energy relaxation time of Tj w 590 ns and 
a dephasing time of about T 2 ~ 500 ns. The sample 
temperature was kept at 33 mK. 


APPENDIX B: ADDITIONAL EXPERIMENTAL 
DATA 

Figure Al shows two examples of measured relaxation 
curves of the 3D-transmon qubit and the fits to the data. 
We fit the measurements to decay curves of the form 
Ae~ Tlt + B with the free parameters A, B and Ti. We 
show one trace where the fit converged with a very small 
standard error (a) and another where the convergence 
was worse (b). The second trace might be better fit by 
assuming a double exponential decay where at some time 
the decay rate changed spontaneously due to a change in 
the environmental noise spectrum (not shown). 

Figs. A2-A4 show the full datasets of the fluctuations 
in the relaxation rate T i measured in our 3D-transmon at 
three different experimental temperatures. We also show 
the histograms for the probability of occurrence of a par¬ 
ticular value of Ti for all three temperatures as well as 
the fluctuations in the fit amplitude A and background B. 
The later two show some fluctuations, but are relatively 



time [ms] 

Figure Al. (Color online) Examples of decay curves taken at 
30 mK. Points are experimental data of the qubit excitation 
probability P(|l)) as a function of time after an initial n- 
pulse applied at t = 0, normalised to lie between 0 and 1. 
The red curve is the result of a fit of the data to the function 
P(|l}) = Ae -rit + B, with free parameters A , B and Ti. We 
show one curve with minimum standard error in the decay 
amplitude A — B (a) and another curve with maximum error 
(b). The lower curve might be better described not by purely 
exponential decay, if for example during measurement of the 
data the noise spectrum shows a sudden jump. 


flat on the scale of the changes observed in Ti. Ampli¬ 
tude fluctuations might be explained if the qubit’s level 
splitting varies in time, which, together with a strongly 
coloured high-frequency noise spectrum provides an al¬ 
ternative model for the fluctuations in the qubit’s relax¬ 
ation rate (c.f. main text). From the data in Figs. A2- 
A4, we conclude that this mechanism might be present 
but is weak and not the main contribution. Additionally, 
we show the two-time correlation function of the relax¬ 
ation rate as well as its Fourier transform. We fit the Tj- 
fluctuation spectrum to two different functions and show 
the results in the plots. The red dashed lines are from 
the best fit to the function d/w“, corresponding to a 1//- 
type frequency distribution as it is expected from a dense 
distribution of low-frequency TLS 9,15 The blue dashed 
lines are results from a fit to a zero-frequency Lorentzian 
~ A'y/ ( 7 2 +w 2 ), as it would result from a single dominant 
low-frequency TLS, c.f. Eq. (2). For our data presented 
here, the temperature dependence of the fluctuation am¬ 
plitude is inconclusive and does not give any indication if 
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Figure A2. (Color online) Experimental data on T\ fluctua¬ 
tions in the 3D-transmon sample at a temperature of 30 mK. 
(a) shows the relaxation rates F i from fits of the experiments 
to an exponential decay curve, P(| 1)) = Ae~ Vlt + B, with er¬ 
ror bars corresponding to the 95% confidence interval of the 
fits. The black dashed lines are a moving average over 10 
points and the red dotted lines are the mean values over the 
full dataset. The inset shows a histogram of the probabili¬ 
ties of values for the relaxation rate T i. (b) shows the time 
evolution of excitation amplitudes A and background B from 
the same fits, including error bars and moving averages in 
black, (c) depicts the absolute value of the Fourier transform 
of the two-time correlation function of the relaxation rates 
(Fi(t)ri(O)} in a log-log-plot, with the inset showing the cor¬ 
relation function itself. The red (blue) dashed curve is the 
result of a fit of the data to a A/a)“-spectrum [Lorentzian 
spectrum A-y/{^ 2 + tu 2 )] with fit parameters A = 0.097 and 
a = 0.58 [A = 0.18 and 7 = 0.34 mHz], for details see text. 


our model is accurate. On the other hand, the frequency 
dependence of the correlations seems to follow roughly a 
1/uj dependence, which can be explained in the terms of 
our model. 



0-0 

0 100 200 300 400 500 600 



w/2tt [mHz] 


Figure A3. (Color online) Same as Fig. A2, data taken from 
experiments performed at a temperature of 50mK. Fit pa¬ 
rameters in (c) are A = 0.079 and a = 0.79 [A = 0.087 and 
7 = 0.052 mHz] for red (blue) dashed line. 


APPENDIX C: MODEL CALCULATIONS 


Here we give additional details on the calculations of 
the mean value and spectrum of the Ti-fluctuations in a 
superconducting circuit due to interactions within a bath 
of spurious background TLS. 
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Figure A4. (Color online) Same as Fig. A2, for a sample 
temperature of lOOinK. Fit parameters in (c) are A = 0.067 
and a = 0.45 [A = 0.36 and 7 = 4.91 mHz] for red (blue) 
dashed line. 


TLS parameter distribution 

Rewriting Eq. (8) in terms of the TLS level-splitting 
E and the mixing angle 6 = arctan A/e, we find 

pos 01 0 

P(E, 9)dEd9 = AE a - dEdd . (A.l) 

sin# 

When describing the full distribution of TLS for all ener¬ 
gies, we integrate the tunnel splitting A between AMin it) 

0 and Am 3 x and the asymmetry energy e between eMin = 

0 and eMax- We find for the integration bounds in the 
new variables: # M in = arctan A M in/£Ma x > 0, ^Max = 
arctanA Ma x/£Min = 7t/2 and F M in = \J A^ in + e^ in = 

A M in, -^Max = x/A^ax + e Max-. Here : A M in is defined 
by the minimum tunneling barrier below which the de¬ 
scription as a two-level system breaks down and i^Max 
provides an upper bound on the TLS level-splitting. 

As an example for the distribution of the inter-TLS 
coupling strength g , we write the probability distribution 
in the case where the interaction is mediated by dipolar 
interaction with \g\ ~ 1/r 3 . One finds 

dr _ 4 

P(g)dg = P(r)—dg = p 0 \g\ 3 dg , (A.2) 

og 

where we assumed a constant TLS density in space p 0 . 

Calculating the average 

We here give some of the intermediate steps of the 
calculations of the average qubit relaxation rate and rate 
fluctuations spectrum. 

With the thermal occupation of a TLS in equilibrium, 
(a z ) = cos 6 (a z ) = cos#tanh (E/2kBT), we can directly 
write the mean value of the qubit relaxation rate to lowest 
order in the inter-TLS coupling strength g as 

{%) = 7g 0) + 7^ 5Z 9j cos 6j tanh 

3 

= 7 q 0) + "/q 1] J dg do dE P{g , #, E)g cos 9 tanh , 

(AA) 

where the sum includes all other two-level defects that a 
single high-frequency TS is interacting with. Due to the 
symmetric distribution in inter-TLS coupling strength g , 
the above integral will evaluate to zero. 

For the spectrum of fluctuations of the relaxation rate 
we then calculate the Fourier transform of the rate cor¬ 
relation function as 
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( 7 ^) J29j cos2 °J 


1 — tanh 


2 1 Ej 


2 7i, 


2k B Tj\ 7 i,j 2 + < 


JdgMdEdrnPME.-ntfcofe 


1 - (air) 


271 


7i 2 + w 2 


(A.4) 


where, in evaluating the correlator, we restrict ourselves 
to only the low frequency contribution of the TLS auto¬ 
correlation function Eq. (2), i.e. we focus on TFs with 
E <C k B T. Additionally we are assuming that different 
TLS are uncorrelated, {<J z> j<J z ,i) = 0 and, since we are 
only interested in the bare fluctuations of the rate, we 
have already subtracted the mean rate. 

Finally, adopting the standard assumption for tun¬ 
nelling TLS, P( 71 ) ~ I/ 71 , the frequency dependence 
of the fluctuation spectrum is determined by 



c?7i P(7i) 


271 

7 2 + w 2 


2 arctan 

_ CO 

UJ 


oc 



^ ^ TMax 
^ ^ TMax 


(A.5) 


Here the maximum relevant switching rate 7 Max is given 
by the time of a single Tf-measurement. All faster 
fluctuations will be averaged out in the observations, 
leading to the behaviour ~ 1 /w for u < 7 Max- In the 
opposite case w 7 Max> i-e. when we observe the 
fluctuations on time scales that are short compared to 
l/ 7 Max: the spectrum will show a 1 /w 2 dependence. 


Phonon induced TS switching 


Alternatively to the generic thermally activated 
switching mechanism discussed previously, one can as¬ 
sume a microscopic model for the TLS relaxation rate 
71 . For example for coupling to phonons, and omitting 
irrelevant prefactors, one arrives at' 


71 oc A 2 Pcoth 


\2k B T J 


oc 2 TE 2 sin 2 9 , 


(A. 6 ) 


where in the second step we already assumed that the 
relevant energies of the switching TF are smaller than 
temperature, E <C T. Since the relaxation rate in this 
expression depends mainly on the TF mixing angle 9 , 
the restriction on small switching rates will be realised 
by confining 9 to small values around zero, effectively 
limiting the value of the coupling strength between the 


relevant low-frequency TF and their phonon bath. Phys¬ 
ically, Eq. (A. 6 ) implies that phonons do not induce 
switching in TLS with small tunnelling matrix element 
A. Performing the energy integration in the average we 
then get 


4 TE 2 sin 2 9 
4T 2 E A sin 4 9 + u 2 
4T “+ 4 sin 2 9 
(3 + a)ui 2 ’ 



(A.7) 


where we expanded the integral to leading power in tem¬ 
perature T. Combined with the previous results for the 
prefactor jq 1 \ this leads to the overall temperature and 
frequency dependence of the relaxation rate correlator 


(7,,i(0)7 ? ,iW) u ocw 2 


( rp-^a-2 


—3a 


T 

rj-i 2a+5 


Sto < 72 ,i 
Su ~ 72,i • 


Su> » 72 ,i 

(A. 8 ) 


For the remaining integration over the mixing angle, 
one finds 



d9P{9) sin 2 9 



d9 sin 9 cos“ 9 


1 - COS 1+a 9 M ax 
1 H - OL 


(A.9) 


where $Max is determined from Eq. (A. 6 ) and the value 
of the maximum observable switching rate 7 Max- In 
this case the temperature and frequency dependence to 
leading order in temperature is entirely contained in 
Eq. (A.7). 


Effective inter-TLS interaction range 

Here we give a rough estimate of the maximum inter- 
TLS distance which still allows noticeable interactions 
between them. We assume the TLS to be realised as 
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microscopic electric dipoles of uniform dipole size di = 
le x 10~ 10 m, where e is the charge of a single elec¬ 
tron 5,18 . Then, assuming parallel orientation of the two 
TLS and using the relation between dipole magnitude 
and coupling strength 

2 = 47 r£o£r ( dl, - LC?2, - L “ 3d i,|| rf 2,||) , (A.10) 

we can estimate the maximum distance to effect a mini¬ 
mum coupling strength of c/Min = 1 MHz (c.f. Fig. 3) as 
fMax ~ 110 x 1(T 9 m. The volume in which TLS are in¬ 
teracting strongly enough is thus Vtls ~ 5.6 x 10 -21 m 3 . 
Assuming an overall TLS density of 10 2 /(/xm 3 GHz) 6,13 , 
this leads to the effective frequency density of TLS in the 
interaction region of a single TLS of p ~ 10 _ 1 /GHz. We 
note that the density obtained in Refs. 6 and 13 refers 
only to high-frequency TS, and a much higher density is 
expected for low-frequency TF 5,15 . 


APPENDIX D: RABI-SPECTROSCOPY 


rotating frame at the drive frequency, we then find the 
decoherence rates as 


T-i • 2 o ,1 2 a Sx±_ iy^d) 

r v = sin /3 + - cos /? 71 -—-—- , 

2 SxAUq) 

x 1 2 a 1 1 n ■ a \2 Sx ± (wd + tt) 

r 4 =oCos /? 7n + t (1 — sin/3) 71 ——-—— 
^ 4 bx ± {u q ) 

r 1 2 a , 1 . ■ p\2 S X± (uJ d - H) 

r t=o cos + -(l+sm(3) 71 ——-—— 

£ 4 b X± {U q ) 


(A.13) 


where we defined the rates 

'y v > = \bpx I, (°) , 

71 = ^b 2 ± S x± (^q)y 


7fi = 2 b \\ Sx II ( fi ) 


(A.14) 


and we used the Rabi-frequency f l = + (ui q — uj d ) 2 . 

Here, we introduced the symmetrized correlation func¬ 
tions for the bath variables X, defined as 


When using non-frequency-tuneable qubits like single 
junction transmons, it is still possible to probe the form of 
the noise spectrum in close spectral vicinity of the qubit 
transition frequency. To this end one can make use of the 
fact that for a driven system, the frequencies of the noise 
spectrum relevant for decoherence will be shifted by the 
applied driving frequency. This effect can be thought of 
as a result of interaction of the dressing of the system 
states with drive photons, or similarly in the context of 
sideband transitions. The following derivation is based 
on the work in Ref. 39 , more details can be found there. 

For a two-level system driven with Rabi driving 
strength Ho at frequency u> d we write the Hamiltonian 

H = i uj q a z + H 0 cos u d t a x + Hsys-B + H B , (A.11) 

with the qubit level-splitting u q , the bare Rabi frequency 
fi 0 and driving frequency ui d . For the system-bath cou¬ 
pling term, we take 

Hsy S -B = -b\\(7 z X\\ + -b±a x X± , (A.12) 

where the qubit level splitting is coupled to the bath vari¬ 
able Xy with coupling strength fey and additionally the 
bath variable X± might induce transitions between the 
qubit states due to its coupling with strength b±. Here 
the bath coupling constants b are assumed to be small 
with respect to the other energies in the problem, such 
that we can use perturbation theory in the strength of the 
system-bath coupling term Hsys-B■ We will not specify 
the exact form of the bath Hamiltonian H B but simply 
assume that is of a suitable form to induce Markovian 
decoherence, i.e. it possesses a very large number of de¬ 
grees of freedom and equilibrates on a time scale that is 
much shorter than all system time scales. Moving into a 


S x (u) = ^ (C x (co) + C x {-u)) . (A.15) 

where Cx{u) = J^ x dt e~' UT (X(t)X( 0)) th and the av¬ 
erage (..) th is over the steady state of the bath. For an 
environment in thermal equilibrium, the unsymmetrized 
noise spectrum will follow a detailed balance relation, 
C x (oj) = e~^ ul C x (uj), with the inverse temperature 
/? = 1/ksT. The angle (3 in these expressions defines 
the relationship between drive strength Ho and detun¬ 
ing between drive frequency and qubit splitting and is 
defined as tan/? = Ho/(w g —<*;<*). 

The two rates 7 V and 71 can be determined in inde¬ 
pendent experiments, measuring relaxation from decay 
of the qubit excited state and decay of Ramsey fringes. 
70 on the other hand can potentially be estimated us¬ 
ing 7 ^ and assuming a l//-type dependence of the low- 
frequency noise spectrum. 

In a Rabi experiment, the decay of the oscillations will 
be proportional to e _r2 * with F 2 = F^ + i(F| + Fy.) and 
thus measurements of the Rabi oscillations at different 
drive strengths can be used to infer the noise spectrum 
in the vicinity of the qubit transition frequency u> q . 


APPENDIX E: QUASIPARTICLE DENSITY 
FLUCTUATIONS 

Experimentally it was found that the temperature 
dependence of the relaxation rates of superconducting 
qubits could be well explained when assuming interacting 
with thermally excited quasiparticles 40 . In this theory, 
the low temperature limit of the relaxation time Tj stems 
from assuming a remaining density of non-equilibrium 
quasiparticles, the origin of which is not yet understood. 
Following the ideas developed in Ref. 40 , we conjecture 



14 


that a fluctuating quasiparticle density, i.e. due to re¬ 
combination events or tunnelling to an outside reservoir, 
might lead to the observed fluctuations in relaxation time 
T\. We calculate the required fractional changes in den¬ 
sity as well as in terms of absolute number of quasipar¬ 
ticles for a given qubit design. 

The following calculations follow closely the theory of 
Ref. 40 , and we here only repeat their main steps for clar¬ 
ity. To derive the effects of the interaction between quasi¬ 
particles and superconducting circuits, we start with a 
low-energy Hamiltonian describing tunnelling of quasi¬ 
particles across a Josephson junction at phase difference 
V 


H T = i t ^ sin | a n,J a m,a + h - c - ( A - 16 ) 

n,m,a 


where t is the tunnelling amplitude and the operators 
(inY destroy a quasiparticle in state n with spin cr in the 
left/right lead. Eq. (A. 16) is valid as long as the qubit 
energy u as well as the characteristic energy 6E of the 
quasiparticles is much smaller than the superconducting 
gap A sc , a condition which is well satisfied in experi¬ 
ments. Starting from this equation, the authors in Ref. 40 
derive the quasiparticle linear response function and thus 
the complex admittance of the Josephson junction due to 
quasiparticle tunnelling. 

Using the golden rule, we write the transition rates 
between qubit states due to quasiparticle tunnelling as 


r 


*-»/ — 


(*|sin| |/) 


2 


(A.17) 


Cooper pairs x qp = n qp / 2t'oA sc . Here vq is the density 
of states of electrons in the leads, which we assume to be 
the same on both sides. 

For the case relevant to experiments, where a single 
junction 3D-transmon was used, the relaxation rate due 
to quasiparticles can then be calculated as 


X qp 

1 1-vO - - ~Z~ 

2ir 



(A.19) 


with the junction’s plasma frequency co p = \/8EjEc and 
its charging energy Eq- Eq. (A.19) directly relates a 
qubit’s relaxation rate to the density of quasiparticles. 
From Eq. (A.19), we can extract the fractional quasipar¬ 
ticle density x qp , with the value for the superconducting 
gap of thin-film aluminium 14 : 

A/2t r « 200 /xeV « 50 GHz « 3.2 x 10” 23 J . (A.20) 


Then, for a relaxation time of Tj = 100/us, corresponding 
to Ti = 10 x 10 3 /s (c.f. Fig. 1 in the main text) we find 
the canonical value of x qp «5x 10“ 7 . 

We want to use the relative quasiparticle density de¬ 
termined above to calculate the actual number of quasi¬ 
particles interacting with the qubit sample. For this 
we need the electron density of states at the Fermi 
edge for aluminum, which we take from literature as 
vq = 4.65 x 10 47 m -3 J -147 . We thus find the quasipar¬ 
ticle volume density for the above used relaxation rate 
Ti = 10 x 10 3 /s as 


where u>if = oji — uif is the energy splitting between qubit 
states |i) and |/} and S qp (ui) is the quasiparticle spectral 
density, which can be calculated from the complex admit¬ 
tance via the fluctuation-dissipation theorem. For low 
temperature, T <C A sc and high frequencies 6E, 

one finds 

c t \~ / 2 A sc 

b qp (u)) « x qp - \ - (A.18) 

7T V CO 


n qp = 2v 0 A sc x qp «5x 10 18 m 3 . ( A -21) 

The 3D-transmon used in the experiments consists of 
two paddles of dimensions 350 x 10 -6 • 700 x 10~ 6 • 120 x 
10~ 9 m 3 , with a total volume of Vai ~ 3 x 10” 14 m 3 . We 
then find 

6N qp = V A i n qp /10 w 1.5 x 10 4 (A.22) 


with the junction’s Josephson energy Ej and the frac- as the number of quasiparticles that, for the sample used, 
tional quasiparticle density normalized to the density of leads to a change in the relaxation rate of JT! = lx 10 3 /s. 



